#include <iostream>
#include <fstream>

int main(){
    int n = 128;
    std::ofstream out1("./test/3.data");
    for(int i = 0;i < 4; i = i + 1){
        for(int j = 0;j < 2; j = j + 1){  
            for(int k = 0;k < 2;k = k +1){
                for(int t = 0;t < 2;t = t + 1){
                    out1 << n << "\n";
                    if(j) out1 << "full-weighting" << "\n";
                    else out1 << "injection" << "\n";
                    if(k) out1 << "linear" << "\n";
                    else out1 << "quadratic" << "\n";
                    if(t){
                        if(i==3&&j==1&&k==1)
                            out1 << "V-cycle";
                        else out1 << "V-cycle" << "\n";
                    }
                    else out1 << "full multigrid cycle" << "\n";
                }
            }
        }
        n = n * 2;
    }
    out1.close();
    n = 128;
    double eps = 2.2e-10;
    std::ofstream out2("./test/4.data");
    for(int i = 0;i < 4;i = i + 1){
        for(int j = 0;j < 4;j = j + 1){
            out2 << n << "\n";
            if(i==3&&j==3)
                out2 << eps;
            else out2 << eps << "\n";
            eps = eps / 100;
        }
        eps = 2.2e-10;
        n = n * 2;
    }
    out2.close();
}
